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Quantum chaos in QCD at finite temperature 
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We study complete eigenvalue spectra of the staggered Dirac matrix in quenched QCD on a 6 J x 4 lattice. 
In particular, we investigate the nearest-neighbor spacing distribution P(s) for various values of j3 both in the 
confinement and deconfinement phase. In both phases except far into the deconfinement region, the data agree 
with the Wigner surmise of random matrix theory which is indicative of quantum chaos. No signs of a transition 
to Poisson regularity are found, and the reasons for this result are discussed. 



1. Introduction 

The fluctuation properties of the eigenvalues of 
the lattice Dirac operator have attracted much 
attention recently. In Ref. Q , it was first demon- 
strated that they are described by random ma- 
trix theory (RMT). In particular, the so-called 
nearest- neighbor spacing distribution P(s), i.e., 
the distribution of spacings s between adjacent 
eigenvalues, agrees with the Wigner surmise of 
RMT. According to a conjecture by Bohigas 
quantum systems whose classical counterparts are 
chaotic have a P(s) given by RMT whereas sys- 
tems whose classical counterparts are integrable 
obey a Poisson distribution, P(s) = e~ s . In this 
sense, the form of P(s) characterizes quantum 
chaos. 

At the localization transition in disordered 
mesoscopic systems, one observes a transition in 
P(s) from Wigner to Poisson behavior. The ques- 
tion of whether there is such a transition in the 
case of the lattice Dirac operator was first raised 
in Ref. (!]]. The present study serves to investi- 
gate this question. Recently, a Wigner to Poisson 
transition was also studied in the context of a spa- 
tially homogeneous Yang-Mills-Higgs system j| . 

The lattice Dirac operator falls into one of three 
symmetry classes corresponding to the three chi- 
ral Gaussian ensembles of RMT Q. In Ref. 0, 
the Dirac matrix was studied for SU(2) using 
both staggered and Wilson fermions which corre- 



spond to the chiral Gaussian symplectic (chGSE) 
and orthogonal (chGOE) ensemble, respectively. 
Here, we study SU(3) which for both staggered 
and Wilson fermions corresponds to the chiral 
Gaussian unitary ensemble (chGUE). We thus 
cover the last remaining symmetry class. The 
RMT result for P(s) is quite complicated; it 
can be expressed in terms of so-called prolate 
spheroidal functions, see Ref. Q where P(s) has 
also been tabulated. A very good approximation 
to P(s) is provided by 

P« = ^e~^ (1) 
which is the Wigner surmise for the chGUE. 

2. Eigenvalue analysis 

We generated gauge field configurations using 
the standard Wilson plaquette action for SU(3) 
and constructed the matrix of the Dirac opera- 
tor using the Kogut-Susskind prescription. The 
Dirac matrix is anti-hermitian so that all eigenval- 
ues are imaginary and occur in pairs with oppo- 
site sign. We have worked on a 6 3 x 4 lattice with 
various values of (3 and typically produced 10 in- 
dependent configurations for each f3. All spectra 
were checked against the analytical sum rules 



A„ = 



and 



E A « = 3V > 



(2) 
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where V is the lattice volume. 
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To construct P(s), one first has to "unfold" the 
spectra . This procedure is a local rescaling of 
the energy scale so that the mean level spacing is 
equal to unity on the unfolded scale. Ensemble 
and spectral averages (the latter is possible be- 
cause of the spectral ergodicity property of RMT) 
are only meaningful after unfolding. 

3. Results 

Our results are displayed in the plots. Figure |l| 
shows the staircase function N(E) which is de- 
fined as the number of eigenvalues with A < E. In 
the continuum, N(E) = J Q p(X)d\, where p(A) 
is the spectral density. Note that p(A) cannot be 
obtained from a random-matrix model. 

Figure || shows the nearest-neighbor spacing 
distribution P(s) compared with the RMT re- 
sult. In the confinement phase, we find the ex- 
pected agreement of P(s) with the Wigner sur- 
mise of Eq. ([|) . In the deconfinement phase, we 
still observe agreement with the RMT result up 
to j3 = 10.0. No signs for a transition to Pois- 
son behavior are found. Thus, the deconfinement 
phase transition does not seem to coincide with a 
transition in the spacing distribution. For larger 
values of j3 the eigenvalues start to approach the 
degenerate eigenvalues of the free theory, given by 
A 2 = Yju=i sin. 2 (27rnj,/.Z^,)/a 2 , where a is the lat- 
tice constant and = 0, . . . , — 1. In this case, 
the nearest-neighbor spacing distribution is triv- 
ial. However, it is possible to lift the degeneracies 
of the free eigenvalues using an asymmetric lat- 
tice where L x , L y , etc. are relative primes. For 
large lattices, the nearest-neighbor spacing dis- 
tribution of the non-degenerate free eigenvalues 
is then given by the Poisson distribution. While 
it may be interesting to search for a Wigner to 
Poisson transition on such asymmetric lattices, it 
seems clear that this transition will not coincide 
with the deconfinement phase transition. 

Furthermore, we do not think that the ab- 
sence of a signature for a transition from Wigner 
to Poisson behavior at the deconfinement phase 
transition is due to the finite lattice size. Even for 
the small lattice size we used, the agreement of 
P(s) with the RMT curve is nearly perfect. This 
leads us to believe that we should have seen some 



sign of a transition if it existed in the thermody- 
namic limit. 

4. Conclusions 

We have searched for a transition in the near- 
est-neighbor spacing distribution P(s) from Wig- 
ner to Poisson behavior across the deconfinement 
phase transition. Such a transition exists, e.g., 
at the localization transition in disordered meso- 
scopic systems. In a Yang-Mills-Higgs system a 
smooth transition along a Brody distribution was 
seen || . We found no signature of a transition in 
our lattice data. The data agree with the RMT 
result in both the confinement and the deconfine- 
ment phase except for extremely large values of 
P where the eigenvalues are trivial. 

In hindsight, these results are not totally unex- 
pected. Temporal monopole currents survive the 
deconfinement phase transition leading to con- 
finement of spatial Wilson loops. Thus, even in 
the deconfinement phase, the gauge fields retain 
a certain degree of randomness. In future investi- 
gations one might try to disentangle those spatial 
contributions to the Dirac matrix. 
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Confinement: (3 = 2.8 (solid line) 
and /3 = 5.0 (dashed line) 
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Deconfinement: (3 — 6.0 (solid line) 
and (3 ~ 10.0 (dashed line) 
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Figure 1. The staircase function N(E) representing the number of positive eigenvalues < E for typical 
configurations on a 6 3 x 4 lattice. 
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Figure 2. The nearest-neighbor spacing distribution P(s) averaged over 10 independent configurations 
on a 6 3 x 4 lattice (histograms) compared with the random-matrix result (solid lines). For comparison, 
the Poisson distribution P(s) = e~ s is inserted for f3 — 10.0 (dashed line). There are no changes in P(s) 
across the deconfinement phase transition. 



